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■ Integrals for the product of unitary-matrix elements over the U (n) group 

will be discussed. A group-theoretical formula is available to convert them into 
. a multiple sum, but unfortunately the sums are often tedious to compute. In 

this paper, we develop an alternative method in which these sums are avoided, 
and group theory is rendered unnecessary. Only unitarity and the invariance 



> 

(N 
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. of the Haar measure are required for the computation. The method can also 

o . 

. be used to get a closed expression for the simpler integral of monomials over 

^ . a hypersphere. 

■ I. INTRODUCTION 



X 



The integral 



for the product of n x n unitary matrix elements and their generating functions are useful 
in many areas of physics. That includes two-dimensional quantum gravity [1], QCD, matrix 
models, and statistical and condensed-matter problems of various sorts [2]. These integrals 
are also useful in the parton saturation problem at small Feynman-x [3]. The monomial 
integral above can be computed using a graphical method [4]. A more powerful expression 
can be derived using the Itzykson-Zuber formula [5] as a generating function, or directly 
from group theory [6] using the Frobenius formula [7]. Simplification can be obtained for 
n ^ oo [8], but for finite n, the expression is quite complicated (see eq. (Al)). It involves 
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multiple sums over an expression containing characters of the symmetric group Sp, as well 
as the dimensions of irreducible representations of Sp and of the unitary group U (n) . One of 
the sums is taken over all the relevant irreducible representations, and the others are taken 
over the symmetry groups of the index sets / = {h- ■ - ip) and J = {ji ■ • ■ jp). These sums 
could be long and tedious for large p, and for indices which have a high degree of symmetry. 

In this paper, we discuss an alternative method to calculate the monomial integral, using 
as input only the unitarity of the matrices U, and the invariance of the Haar measure dU. 
No knowledge of group theory is required, and these complicated multiple sums are avoided. 
We shall refer to this method as the 'invariant method'. 

Invariance of the group measure produces various relations between the different in- 
tegrals, which will be discussed in Sec. Ill A. The values of these integrals are obtained 
recursively from the unitarity relation, and that is discussed in Sees. Ill B to III D. 

The invariant method is also applicable to the much simpler case of a monomial integral 
over a hypersphere. This simpler case will be used as a testing ground for the idea. It will be 
discussed in Sec. II, as a preparation for the computation of the unitary integral in Sec. III. 

In a forthcoming paper, the relative advantages of the invariant method and the group- 
theretical formula will be discussed. We will also examine relations that can be obtained by 
combining both approaches. 

II. INTEGRATION OVER A HYPERSPHERE 

Let Jln-i be the unit sphere in n dimensions, defined by 

n 

E^Ni, (1) 

i=l 

and {du)) be its rotationally symmetric volume element, normalized to J{duj) — 1. We 
wish to calculate the integral J{duj)Y over Vtn-i-i where F is a monomial of the coordinates 
Xj (j = 1, 2, • • • , n). This integral is zero unless the power of every Xj is even, in which case 
it can be written in the form 
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(J|J)= j{dw)XjXj, (2) 

where Xj = Y[a=i ^ja is a monomial of degree p, indexed by the set J = (jij2 • • • jp) . 

One can attempt to calculate the integral in several ways. Three of them having analogs 
with the U {n) integrals will be singled out, because the simpler setting of a sphere should 
make their relative merits more transparent. The first two are standard, both using the 
spherical coordinates to calculate. The third one, which we wish to develop in this paper, 
is an invariant approach, requiring no coordinate system in its computation. 

A. Spherical coordinates 

The spherical coordinates of Xj on the unit sphere are: 

Xn = COsOi 

= sin^i cos ^2 

x„_2 = sin 6*1 sin 6*2 cos 6*3 



/n-2 \ 

■^2 = I H sin 9i J cos (j) 

/n-2 \ 

= I sin 9i I sin 0. (3) 

The range of 9i is between and tt, and the range of 4> is between and 27r. The volume 
element is 

{duj') = (^{sin9i)^-'-^d9}j d(j) 



Using the formula 



Jo 



d9 (sin 9y-\cos 9y-^ = - V' ^ / (5) 



2 r(i(r + s 
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and (3), the integral (J| J) can be calculated for every index set J. 

For example, if all the indices in J are equal to n, i.e., J = {nnn ■ ■ - n) = (n^), then the 
integral is equal to 



However, if we replace J = (n^) by J — (1^), namely, replacing in the integrand by 
x'l', then the integral becomes much more harder to calculate, because (n — 1) times more 
integrations must be performed. Yet, on account of the spherical symmetry, the result must 
come out to be the same as (6). This complication arises because a choice of axes breaks 
the spherical symmetry. It can be avoided in the invariant approach discussed below. 

This method relies on an explicit parametrization of fln-i via the spherical coordinates, 
as well as formula (5) allowing the integrations to be carried out. Both become much more 
difficult in the U (n) case, so much so that this method is really not very useful there. For 
that reason there shall be no futher discussion of this method. 



Alternatively, Xj can be expanded in terms of spherical harmonics, and the integral can 
be transformed into a sum using the orthonormality of the spherical harmonics. For n = 3, 
the expansion is 



{J\J)^S{p) 




/o^/'(cosgi)^P(singi)^-^dgi 
/o"/'(sin^i)"-2d^i 



(6) 



B. Group theory 



(7) 



and the integral becomes 



(8) 



For n > 3, many more sums are involved in eq. (8). 
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There are two non-trivial tasks in this approach: to find the coefficients aem, and to carry 
out the sum in (8). These tasks become quite difficult in practice for large n or p. 

There is an analogous group-theoretical technique to calculate the U (n) integral, which 
is reviewed in Appendix A. Using the Probenius relation, or the Itzykson-Zuber formula, 
a formula can be derived to turn the integral into a multiple sum. As mentioned in the 
Introduction, the sums could be very involved, so in practice this method may not be the 
best way to obtain a result. The invariant approach discussed below and in the next section 
might be simpler. 

C. The invariant approach 

The integral in (2) can be calculated directly, using only condition (1) and the invariance 
of (du) under rotation. In particular, there is no need to employ the spherical coordinate 
system, and no need to know any integration formula, nor group theory. 

The invariant approach will be used in the next section to calculate integrals (of mono- 
mials of unitary matrix elements) over the unitary group U{n). In that case, (1) is replaced 
by the unitarity condition, and dcu is replaced by the invariant Haar measure dU of the 
unitary group. 

It is convenient to arrange the p indices in J according to the distinct values (between 
1 and n) they take. If mi of these p indices take on a value, m2 of them take on a second 
value, and so on, then the integral {J\J) will be denoted by S{mim2 • ■ -mt), where t is the 
number of non-vanishing mj's, and J2l=i rui — p. Spherical symmetry guarantees that the 
integral is independent of the the specific values the indices assume. This means, among 
other things, that S is symmetrical in all its arguments. 

To calculate S, the invariance of du is used to relate the various 5"s to S{p). Then the 
value of S{p) is calculated using the sphere condition (1). 

The invariance of duj can be exploited in the following way. A rotation in the XiXk plane 
by an angle ^, 
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Xk -sxi + cxk, (9) 

where c = cos^ and s = sin^, will leave dcu invariant. Equivalently, if we subject the 
integrand XjXj in ( J| J) to such a rotation, the integral S'(mi • • • ruk) will remain unchanged. 

Let us start out with the integral S{p) whose Xj is equal to xf. Under (9), with (i, k) — 
(1, 2), the integrand becomes 

{xlY ^ (cxi + sx2f^ = E (^2^) {cx,f^'-'\sx2f' + ■■■. (10) 

The ellipsis indicates terms odd in xi and X2, which can be dropped because they do not 
contribute to the integral. The invariance of the integral under this transformation then 
yields the relation 

SiP) = E (2^) c^-'-'^'^'^Sip - e, e). (11) 
Since this is true for all ^, the right-hand-side must be independent of ^. That requires 

S{p-e,e) = j^S{p). (12) 
Similarly, we can apply (9) and the whole procedure to (i, k) — (2, 3) to get 

(e) (p) (e\ 

S{P - e, e - /, /) = ji^Sip - e, e) = T^j&^iP)- (13) 

(S) (X) (5) 

Continuing thus, we finally obtain 

c, , (gS.i ".)!(ni.i(2r».)!) ^, , ,,,, 

To complete the calculation we must calculate S (p) . This can be done by using condition 
(1). Since S{p — 1, 1) is independent of what coordinate xj the multiplicity 1 sits on, as long 
as it is not on the coordinate whose multiplicity is p — 1, the sphere condition (1) can be 
translated to read 

_ i)s{p - 1, 1) + S{p) = S{p - 1). (15) 



Using (14), we know that S{p — 1, 1) = S{p)/{2p — 1). Substituting this back into (15), we 
get the recursion relation 

^W = ;^|^so,-i). (16) 

With the initial value S{0) — 1, (16) can be solved to yield 

(2p-l)(2p-3)---l _ 1 r(p+^)r(in) 
""^^ {n + 2p-2)in + 2p-A)...n ^ v(p^\n) ' ^ ^ 

which agrees with the answer given by (6). 

The general result is obtained by substituting (17) into (14). 



III. INTEGRATION OVER THE UNITARY GROUP 

Let Uij denote the (ij) matrix element of an n x n unitary matrix, and U*j its complex 
conjugate. The product na=i ^iaia abbreviated as C/jj, with the index sets being I — 
• • - ip) and J = {jij2 • • • jp)- We shall refer to p as the degree of Uij. Since the matrix 
elements commute with one another, the order of the indices is irrelevant. Thus if P e S'p 
denotes a permutation of the p indices, and the permuted index set is denoted by Ip — 
{ip{i)ip{2) ■ ■ ■ ip{p)), then 

Uij^Ui.j,. (18) 

We want to calculate the monomial integral 

{I,J\K,L)^ j{dU)U*jUKL (19) 

over the unitary group U{n). The degree of Ujj is assumed to be p and that of Ukl is 
assumed to be q. The Haar measure dU appearing in (19) is left- and right- invariant, and 
normalized to J{dU) — 1. 

As mentioned in the last section, a group-theoretical formula to calculate the integral is 
available (see Appendix A). The integral is expressed as a multiple sum, with a summand 
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involving the character of the symmetric group Sp, the dimensions of the irreducible rep- 
resentations of U (n) and Sp, as well as the index structure of /, J, K, L. This formula is 
general, though not always the best way to obtain the result, because the multiple sums 
are often tedious and difficult to do. In what follows, we shall develop another method to 
calculate the integral, using the invariant approach discussed in the last section. With this 
approach, no knowledge of group theory is required, and multiple sums are avoided. 
The invariant calculation relies only on the unitarity of the matrices in the integrand, 

n n 

E u^Ju;^ = E u,^u;^ = Su, (20) 

i=i i=i 

as well as the invariance of the Haar measure. The latter implies 

J {dU)f{U, U*) = J {dU)f{VU, V*U*) = J {dU)f{UV, U*V*) (21) 
for any function /, and any V e U{ri). We shall apply (21) to the function f{U,U*) — 

The calculation is very similar to that in the last section, though more complicated. The 
spherical condition (1) is now replaced by the unitarity condition (20), and the rotational 
invariance of dcu is now replaced by the group invariance of dU. 



A. Relations from invariance 

Eq. (21) is very powerful. Depending on the choice of V, many relations can be derived. 
Here are some examples. 

1. Vij = e*"^ Sij 

With this choice, f{UV, U*V*) = e'^i-P^^f{U, U*). Hence the integral (19) is zero unless 
p = q. For this reason we shall assume q = p from now on. 
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2. Vi 




With this choice, f{UV, U*V*) = e'^f{U, U*), where ^ = El=Mia - (f>jj- In order for 
the integral (19) not to be zero, this phase ^ must vanish. If all (piS are different, this 
happens only when the indices in L are permutations of the indices in J. Similarly, one can 
show that the indices in K must be permutations of the indices in I. If it! and S are the 
permutations in Sp that do the job, then 



Using (18), we may assume K — I, with the non-zero integrals (19) now being of the form 



for some Q e Sp. 

The non-zero integral over a sphere calculated in the last section is of the form (J| J), 
but the non-zero integral over U{ri) is of the form (7, J|7, Jq). The presence of an additional 
index set I, and the possibility that Jq ^ J, both make it harder to calculate the unitarity 
integral than the spherical integral, though the idea is precisely the same. 

If Jq — J, then the integral (23) is positive definite. We shall refer to integrals of that 
type as direct integrals. If Jq ^ J, the sign is not guaranteed, and we shall refer to those 
integrals as exchange integrals. 

The direct integral looks deceptively similar to the spherical integral (2) in a complex 
n^-dimensional space. By mapping Uij to Xa, with a = running between 1 and n^, the 
complex equivalent of (1), namely Z^^li l^ap = n, is guaranteed by the unitarity relation 
(20). One might therefore think that the direct integrals would turn out to be very similar 
to the spherical integrals, whose result is given by (14) and (17). Unfortunately that is not 
the case, because the measure dU is not rotational invariant in the n^-dimensional complex 
vector space. As a result, even the direct integrals become more difficult to calculate than 
the spherical integrals of the last section. 



K^Ir, L^Js. 



(22) 



{I,J\I,Jq) 



(23) 
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3. permutation 



Now, choose F to be a permutation matrix of n objects. Then VU is obtained from U 
by permuting its rows, and UV is obtained from U by permuting its columns. With this 1/, 
eq. (21) implies that 

{I,J\K,L) = {I',J\K',L) = {I,J'\K,L'), (24) 

where / ' is obtained from / by a reassignment of the values of its indices, and K ' is obtained 
from K by the same reassignment. For instance, let us take p — 6 and n — 8. Suppose 
I = (11555247) and K = (51524571) (recall from (22) that K has to be a permutation 
of /). If we make the reassignment 1 8,2 •e^ 5,4 7, then /' = (88222574) and 
K' = (28257248). 

In other words, the integral is affected by whether the indices take on the same or different 
values, but is independent of what these values are. 

4- row-column interchange 

Since the measure is invariant under transposition, (df/) = {dU^), the integral is unaltered 
if we interchange the rows with the columns: 

{I,J\K,L) = {J,I\L,K). (25) 

5. graphical representation 

It is convenient to employ a graphical description of the integral (/, J\I, Jq), as illustrated 
in Fig. 1. The dots in the left-hand column represent the distinct values in the index set /, 
and the dots in the right-hand column represent the distinct values of the index set J. As 
per (24), the integral does not depend on what these values are. Since Jq is a permutation 
of J, it shares the same distinct values, and hence the same dots as J. 
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Factors of Ujj are shown as (thin) sohd hnes, and factors of Uj^j^ are shown as dotted 
hnes. A number appearing above the sohd hne or below the dotted hne denotes the mul- 
tiphcity {i.e., the power) of that matrix element. If the number is absent, the multiplicity 
is taken to be 1. When a pair U*jUij occurs together, we may choose to replace its pair of 
thin-solid and dotted lines by a thick solid line. In that case the multiplicity designation 
refers to the multiplicity of the pair. Thus a direct integral can always be drawn with only 
thick solid hnes. 




FIG. 1. Graphical representation of the integral {I,J\I,Jq), where / = (1323245), J = (1227^), 
and Jq = (217^27). 

When multiplicity is taken into account, the number of solid lines and the number of 
dotted hnes connected to each dot must be equal. This simply reflects the fact that Jq is a 
permutation of J. 

The integral {I,J\I,Jq) with I = (11123345) = (1323^45), J = (12277777) = (12^7^), 
and Jq = (21777727) = (217^27), is depicted in Fig. 1. The 5 left-hand dots represent the 
five distinct numbers 1,2,3,4,5 in /, and the three dots in the right-hand column represent 
the three distinct numbers 1,2,7 in J and Jq. 

The invariance of (25) means that we may switch the left-hand column of dots with the 
right-hand column of dots. Namely, a reflection about the vertical line halfway between the 
two columns will not change the integral. 
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6. rotation 



Choose V — R{ab) to be the matrix which rotates the (ab) plane by an angle ^. This 
matrix has I's along the main diagonal, except at the (aa) and (bb) positions, where the 
matrix element is c = cos^. The off-diagonal matrix elements are all zero, except at the 
positions (ab) and (ba), where the matrix elements are respectively s = sin^ and —s. 

The replacement U — > UV causes the following change in the matrix elements: 

Uia +C Uia + S Uib, 
Uib -S Uia + C Uib, 

Uij ^ Uij, (26) 

provided j 7^ a,b. Similar replacements on U*^ and [7*^ should also be made. The result is 
to change (/, J\I, Jq) into a sum of terms of the form^ Me(c^)'^~'^(s^)^, where d is the total 
number of column indices in U^j with value a or 6, and e varies between and d. The 
invariance condition (21) requires 

(/,J|/,jQ) = X:Me(cY-(s2)^ (27) 

e=0 

In order for this to be true for all ^, we must have 

Me = (/,J|/,JQ)Q =MoQ, (28) 

where {^^ = d\ / e\{d — e)\ is the binomial coefficient. 

Let us see how Mg is computed in the graphical language. Take any two dots on the 
right-hand column to represent the values a and b. One of the two dots should have some 
lines attached to it, but the other one may or may not be empty. The total number of solid 
(or of dotted) lines attached to the two dots is d. Now move e (thin) solid and e dotted 
hues between the two dots, subject to the constraint that at the end of the move, each dot 



^Odd powers of cs never enters because of (22). 
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must have an equal number of solid and dotted lines attached to it (otherwise the integral 
is zero). Assign a weight +1 for a line moved from a to b, and a weight —1 for each line 
moved from b to a. The quantity Mg in (28) is simply the sum of all the resulting integrals 
after the move, weighted by the product of the factors ±1 associated with each move. 

It is important to note that these relations are local. They involve only the indices of J 
and Jq with values a and b. It does not matter what / is and what the rest of the indices 
of J and Jq are. 

Let us illustrate this graphical application of the rotational relation with two examples, 
example 1 

Fig. 2 represents part of a diagram. The whole diagram may have many more dots and 
lines. They are not drawn because the relation derived below is independent of these other 
dots and lines. 



l2 • 




T« a 




a 



(a) 



(b) 



Is •- 



b 



a 




* b 



a 



I4 



I4 •' 



(c) (d) 
FIG. 2. Four diagrams illustrating how rotational invariance of the measure can be used to relate 
these four integrals. See eq. (29). 
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The values of a and b must be different, and tfiere must be no more lines coming out of 
them than is shown in Fig. 2. The indices ii, ^2, '^3, H may or may not have the same values, 
and there may be many more lines attached to them than is shown, as long as these other 
lines do not connect with either a or b. If some of their values are the same, say ii — i2, 
then graphically the two dots ii,i2 simply merge together into a single dot. If the values 
are different, then other lines must also come out of these dots in order to make the integral 
non-zero. 

The integrals corresponding to Figs. 2(a), 2(b), 2(c), and 2(d) will be denoted by /(2a), 
I {2b), I {2c), and I {2d), respectively. In these diagrams we are dealing with d = 2, because 
there are two pairs of (solid and dotted) lines ending on the right-hand dots. 

Let us apply rotation to Fig. 2(a). We could either move two hues (e = 1), or four lines 
(e = 2). By moving two lines, we get I {2b) (two solid lines), /(2c) (two dotted lines), and 
two others, h (one solid and one dotted hues from a to b) and I2 (one sohd and one dotted 
lines from 6 to a). The graphs for the last two are not shown, but they can be obtained from 
/(2a) by merging the a and b dots. In this way we get Mi = —I{2b) — /(2c) + Ii + I2 = 
2{—I{2b) + Ii), where (24) has been used in the last step. By using (28), we conclude that 
Ml — 2Mo = 2/(2a). Hence we obtain the relation 

/(2a) = -/(26) +/i. (29) 

If we move all four lines, we get M2 — I {2d). The formula in (28) demands M2 — Mq, 
or I {2d) = I {2a). We already know this to be true from (24). 

example 2: the fan relation 

Fig. 3(b) is a partial diagram of some integral. We assume there are no other lines 
connected to the dots on the right-hand column, though there may be other lines emerging 
from dot i. There may also be many other dots and lines not shown in the diagram. 
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(a) f(d) 



(b) f(mj,m2, ... ,mj 



FIG. 3. Using rotational invariance of the measure, the single line in (a) can be spread out into 
the fan of lines in (6) to get the relation displayed in (30) . 

The multiplicities of the lines shown are rrii {1 <i <t), so the integrand is proportional 
to \Uia\'^"^'^\Uib\'^"^^ ■ ■ ■. The corresponding integral is denoted by /(mi, m2, • • • , rrit). 
Using (28) repeatedly, it will be shown below that 

/(mi,m2,---,mt) = ^- r-/(a), (30) 

where d = ELi and f{d) = f{d, 0, 0, • • • , 0). 

f{d) is drawn in Fig. 3(a). Having multiplicity d means that there are d (thin) solid 
lines and d dotted lines between i and a. Now move e = I]*=2^i solid and e dotted lines 
from a to an empty dot b. There are ways of choosing the set of solid lines to move, and 
independently there are also (^fj ways to select the dotted hues. Hence Mg = (^fj f{d — e, e). 
From (28), we know that = QMq = Qf{d). Hence 

f{d-e,e) = -^fid). (31) 
15 



This process can be repeated by moving g — rrii pair of lines from b to an empty 



dot c. Then we get 



f{d- e,e- g,g) 



1 



f{d-e,e) 



1 1 



fid). 



(32) 





By repeating this process again and again, we arrive at the fan relation (30), which tells us 
how to fan out a thick line with a high multiplicity into t different lines. 



In the last subsection, relations between different integrals are obtained using the invari- 
ance requirement (21). To calculate the actual value of any of these integrals, the unitarity 
condition (20) must be used. 

The unitarity sum (20) for i ^ I simply brings out more relations between different 
integrals. But for i — I, a pair of U, U* disappears on the right-hand side of (20), so (20) 
relates integrals of degree p to integrals of degree p — 1. By using this repeatedly, eventually 
the degree comes down to zero, and the integral is known to be 1. In this way the values of 
the integrals can be computed recursively. 

This procedure will be illustrated by various examples in the rest of this section. 



To get this result, we make use of the fact from (24) that (i, a|i, a) is independent of i and 
a. Summing over a (from 1 to n), and using (20), we get n(i,a|i,a) = 1. Hence eq. (33). 
This calculation can be generalized to the integral in Fig. 4(a) to give 



B. Results from unitarity 



1. direct integrals 



the fan integrals 



The simplest {p—1) direct integral is 




(33) 
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IT? — 1 ) ^771 ^ 

((n, inun, in) ^ n^) = ^l^Ji^y - (34) 

Note that Fig. 4(a) is just the diagram Fig. 3(a), but without any additional dots and 
hnes. 



i 




(c) Z(mj,m2,mj) (d) 



FIG. 4. Diagram (6) shows how the integral in (a) can be obtained from unitarity by summing 
over the index indicated by an arrow. A recursion relation whose solution is (34). Similarly, 
unitarity applied to the indicated index in (d) yields a recursion relation whose solution gives (38) 
for the 'Z integral' in (c). 

To obtain this result, start from the integral F{m — 1, 1), shown in Fig. 4(b). The integral 
is independent of the value of the index (dot) indicated by a downward arrow, as long as 
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it does not take on the value of the other dot. Sum over this indicated index, from 1 to n. 
Eq. (20) imphes 

(n- l)F(m- 1,1) + F(m) = F(m- 1). (35) 

Now the fan relation (30) tells us that F(m — 1, 1) = F{m)/m. Substituting this into (35), 
we get a recursion relation between F{m) and F{m — 1), namely, 

rrr 

F(m) = F(m-l) . (36) 

Using the initial value F(0) = 1, this recursion relation can be solved to get (34). 

Define the 'fan integral' F{mi,m2, • • • , m„) to be Fig. 3(b), without any extra dots and 
lines. It follows from (34) and (30) that 

F[mi,m2,---,rnt) = -r r;- (37) 

(n + E-=i rrii-iy. 

the Z integrals 

Next, consider the 'Z integral' in Fig. 4(c). We shall prove that 

Zim m m ) = m^mM^Mn - 2)\{n - l)!(n + mi+m^- 2)! 
^ ^' ^' ^' (n + mi-2)!(n + m3-2)!(n + mi + m2 + m3-l)!' ^ ' 

To do so, consider Fig. 4(d). Summing over the index indicated by the vertical arrow, and 
denoting the integral in Fig. 4(d) by /(4d), the unitarity condition (20) gives 

(n - 2)7(4d) + Z{mi, m2, mz) + Z{mi, 1712 + 1, ma - 1) = Z{mi, m2, ms - 1). (39) 

The fan formula (30) tells us that /(4ci) = Z(mi, 777,2, m3)/m3. Substituting this into (39), 
we get a recursion relation in 7773: 

Z{mi, 7772, "73) = — [Z(7?7i, 7772, 7773 - 1) - Z{mi, 1712 + 1, 77l3 - 1)] . (40) 

n + 7713 ^ 2 

Using the initial value Z{mi, 7772, 0) = F{mi, 7772) — {n — l)!777i 17772!/ (71 + 7771 + 7772 — 1)!, the 
recursion relation can be solved to yield (38). 

Relation (30) can also be used to fan out the two open ends of Fig. 4(c) and (38) to 
obtain the 'fanned Z integrals'. 



C. Exchange integrals 



To illustrate how to compute exchange integrals, all second and third degree exchange 
integrals will be computed in this subsection. 

1. p = 2 

All the second degree integrals are shown in Fig. 5. The integrals in Figs. 5(a), 5(c), 
5(d), 5(e) are direct integrals, either of the fan type, or the Z type, so they are known. The 
only exchange integral is E{2) depicted in Fig. 5(b). It can be computed either by rotation 
from a direct integral, or by unitarity. We will discuss both methods. 



(a) Z( 1,0,1) 




(b) E(2) 





(c) F(l,l) 



(d) F(l,l) 



(e) F(2) 



FIG. 5. Integrals of the second degree. 
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by rotation 

Starting from Fig. 5(d), rotate two solid and two dotted lines from the dot in the right- 
hand column to an empty dot. This is a special case of example 1 of Sec. IIIA5 and Fig. 2, 
but let us do it directly once again. Using the notation in (28), we get 

Ml = 2Z(1, 0, 1) + 2E{2) = 2Mo = 2F(1, 1). (41) 

Hence 

£;(2) = F(1,1)-Z(1,0,1)= , ^ - = , J (42) 

^ ^ ^ ^ V > > y n{n+l) (n-l)(n + l) n{n^ - 1) ^ ' 

by unitarity 

Summing over the indicated index in Fig. 5(b) from 1 to n yields (n— 1)£'(2)-|-F(1, 1) = 0, 
hence 

E{2) = -^^illll = — i -. (43) 

^ ' n-1 n{n^ - 1) ^ ^ 

2. p = 3 

Fig. 6 shows the two direct integrals which are not of the fan type or the Z type, and all 
the exchange integrals of degree 3. 

To get /(6a), the integral for Fig. 6(a), carry out a unitarity sum on the indicated index. 
This yields {n - l)7(6a) + F(l, 1, 1) = F(l, 1), hence 

^^^"^"^ ~ n-1 (n(n + 1) ~ n(n + l)(n + 2)) ~ (n- l)n(n + 2)' ^^^^ 
This integral can also be computed by fanning out the bottom line of Z{2,0, 1) = 2{n + 
l)(n-2)!/(n + 2)!. 

To compute I {6b), take a unitary sum over the indicated index. This yields (n — 2)7(6&)-|- 
2/(6a) = ^(1,0,1). Hence 

^ ^ n-2\n^-l {n-l)n{n + 2)) {n + 2)\ ^ ' 
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FIG. 6. Integrals of the third degree. 



The exchange integrals can be obtained by taking the unitary sum on the indicated 
vertex. Note that the right-hand side of the sum is always zero in the case of exchange 
integrals. 

In this way we obtain the relations 

(n- l)/(6c) + F(1,2) = 
(n- l)/(6ci) + F(1,2) = 
(n- l)/(6e) + F(1,1,1) = 
(n - 2)/(6/) + /(6a) + /(6e) = 
21 



{n - 2)I{6g) + 27(6e) = 0. 



(46) 



The solutions are 



/(6c) = I{6d) 
/(6e) 



2(n- 1)! 



= -2 



(n + 2)!(n-l) (n + 2)! 

F(l,l,l)_ (n-2)! 



n-l 
1 



n-2 
2 

' n-2 



(n + 2)! 

(J(6a) + /(6e)) = 
,(n-3)! 



/(6e) 



(n + 2)! 



(n~3)! 
'(n + 2)! 



(47) 




D. The X integrals 





(a) 



(b) 



(c) 





(d) (e) 

FIG. 7. (a) The 'X integral' and its weights. (6) — (e) Diagrams obtained from (a) by rotating 
away a pair of hnes from the r — s junction. The weights in these four diagrams are the same 
as those in (a), except for the ones with a '— ' sign, in which case the corresponding weight is 
decreased by 1. 
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To illustrate how direct and exchange integrals may be coupled in the recursion relation, 
let us look at the 1-loop 'X integrals' depicted in Fig. 7(a). This integral, specified by 
the four weights r, s, t, u of U* and the four weights r', s', t', u' of t/, will be designated 
as X{rstu\r's't'u') — I {7a). Since the number of dotted lines and the number of solid lines 
emerging from each vertex must be equal, there are three relations for these eight parameters, 

r' + s' — r + s, 
s' + t' ^ s + t, 

t' + u' ^t + u, (48) 

so only five independent parameters are required to specify all the X integrals. 

It is fairly complicated to calculate all these integrals, so we will only derive the recursion 
relation here and illustrate how it can be used in the simplest case. Let p = r + r' be the 
total number of top lines in X, and a — s + s' the total number of lines on one side. The 
idea is to find a recursion relation in p + a, each time reducing either p or cr by 1. Eventually 
one gets down to either p = or cr = 0, which are the Z integrals obtained before. 

1. recursion relation 

Rotate one solid and one dotted line from the r — s junction of Fig. 7(a) to an empty 
dot. The result is Fig. 7(b) to Fig. 7(e). In the notation of (28), we have 

Ml = r'sl(7b) + rs'I(7c) + ss'I(7d) + rr'I(7e) = (r + s)Mo = (r + s)I(7a), (49) 

where I {7b) is the integral depicted in Fig. 7(b), etc. 

The parameters of Figs. 7(b) to 7(e) are those of 7(a), except where a ' — ' sign occurs, 
in which case the corresponding parameter is reduced by 1. In Fig. 8 we will also use a '+' 
sign to indicate where the parameter is increased by 1. 

The unitarity sum, apphed to the indicated index in Figs. 7(b) to 7(e), results in the 
relations 
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= (n - 2)7(76) + I{7a) + 7(86) 
= (n - 2)7(7c) + 7(7a) + 7(8c), 
7(8^2) = (n - 2)7(7d) + 7(7a) + 7(8di), 

/(8e2) = (n - 2)7(7e) + 7(7a) + 7(8ei). (50) 

A substitution of (50) into (49) yields the desired recursion relation 

s r's7(86) + rs'7(8c) + ss'[7(8rfi) - 7(8^2)] + rr'[7(8ei) - 7(862)] , 

liia] — ; . oi 

^ ^ {r + s){r' + s' + n-2) ^ ' 






(b) 



(c) 



(di) 






(dz) (ei) (e2) 

FIG. 8. Terms obtained from the unitarity sum of Figs. 7(6) to 7(e). The weights in these 

diagrams are the same as in Fig. 7(a), except on those hnes indicated by a '— ' or a '+'. In those 

cases, the corresponding weights are decreased or increased by 1. 



Let us illustrate the recursion relation by computing the simplest cases, with p + a — 2. 

There are four possibihties with r + r' + s + s' — 2. Two of them are the known 
Z-integrals, namely, X{0,l,t,u\0,l,t,u) = Z{l,t,u) and X(l, 0, t, -ull, 0, t, m) = Z{l,u,t). 
The other two are exchange integrals which can be obtained from (51): 
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X(1.0.,„|0. M - 1) ^ -iM ^ -^<i^ . - V^!'::^' . (52) 

n — 1 n — 1 [n + t + uy. 

X(0, 1,^-1, 0, - 1) = = = [n + t + u)\ " ^^^^ 

In particular, (52) reduces to /(6c?) of (47) when t — u — 1. 

IV. CONCLUSION 

We have shown how an integral over the U {n) manifold can be computed recursively using 
only the invariance of the Haar measure and the unitarity condition. The same method can 
also be used to get a closed expression for a monomial integral over a unit sphere. 

In a forthcoming paper, we shall compare the relative advantages of the invariant method 
developed here, and the group-theoretical method reviewed in Appendix A. We will show 
how the two methods can be combined to strengthen each other. 

CSL would like to thank G. Semenoff and A. D'Adda for helpful discussions. This 
research is supported by the Natural Sciences and Engineering Research Council of Canada 
and by the Fonds de recherche sur la nature et les technologies of Quebec. 

APPENDIX A: GROUP THEORETICAL CALCULATION 

As shown in (23), the non-zero integrals of (19) can be written in the form (7, J|7, Jq). 
Using group theory, to be reviewed below, the integral can be turned into a multiple sum^: 

(7, J|7, Jq) = E EE T^Xf{SQR\ (Al) 

where p is the degree of Ujj appearing in the integral (19). The other symbols will be 
explained below. 

The irreducible representations of the unitary group U{n) are labeled by a Young's 
tableau. It is defined by a sequence of non- negative integers / = (/1/2 • • • fn), with /j > /j+i. 



^CSL wishes to thank Prof. Alessandro D'Adda for introducing him to this formula. 
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All irreducible representations contained in a pth rank tensor have their Young's tableaux 
restricted to p boxes, namely, ft = p. In that case clearly /, = for i > p. It 
is customary to drop the zeros when the sequence / is written, so / can be written as 
(/1/2 ■ • ■ fp), or even shorter if there are more zeros. 

The dimension of the irreducible representation / of U{n), denoted by df, is given by 
the ratio of two Vandermonde determinants 

= i^(n-l,n-2,...,0)' ^^2) 

where ii = fi + n — i, and 

n 

D{xi,X2r--:Xn) = J| (Xj - ) . (A3) 
i>j=l 

The irreducible representations of the symmetric group Sp are also labeled by Young's 
tableaux / = (/1/2 • • ■ fp)- The dimension of an irreducible representation of Sp is denoted 
hy (if, and the character for the element P & Sp'is denoted by X/(-P). Tables are available 
to give their values for small p. 

The character Xf{P) depends only on the class that P belongs to. If a permutation P 
is written in the cycle form, then permutations with the same cycle structure belong to the 
same class. The cycle structure can be labeled by a Young's tableau c = (ciC2 • • • Cp), where 
Ci is the length of the longest cycle in P, C2 is the length of the next longest cycle in P, 
etc. If P e c, we will also write X/(P) as X/(c). The characters together are given by 7^ 
numbers, where 7p is either the total number of distinct classes in Sp, or the number of 
inequivalent irreducible representations. It is equal to the number of partitions of p, and is 
given by 7^ = 1, 2, 3, 5, 7, 11, 15, for p = 1, 2, 3, 4, 5, 6, 7, respectively. 

Qi C Sp is the symmetry group of the index set /, and Qj C Sp is the symmetry group 
of the index set J. For example, if / = (111338888), then ^ S^ x S2 x S^ C Sg. If 
/ = (13254798), then Qj consists of the identity e only. 

The sum in (Al) is over all the irreducible representations /, all elements R of the 
symmetry group Qi, and all elements S of the symmetry group Qj. 
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The simplest integrals to calculate are those where the indices in I all take on distinct 
values, and similarly for J. In that case, Qi — Gj = e, so the sums in (Al) reduces the 
single sum over the irreducible representations /. Since XfiQ) depends on the class c Q lies 
in, there are 7p distinct integrals (/, J|/, Jq) of this type. It is convenient to denote these 
integrals by ^(c). To compute them, we need to use (A2) to compute d/, a character table 
of Sp to compute Xf ^^'^ df, then we must sum up jp terms in (Al) before we get ^(c). 

More generally, (Al) can be written in terms of ^(c) as 

(7,J|7,jQ)=^Ar(7,J,g|c)e(c), 

c 

where 

N{I,J,Q\c)^ E S{SQRec) (A4) 

ReGi SgGj 

is the total number of SQR in class c. They are often quite tedious to compute. Once it is 
calculated, we still have to carry out the 7^ sums over c to get {I, J\I, Jq). 

This completes the description of formula (Al). In the remainder of this Appendix, we 
will sketch how it is arrived at. 

The orthonormal relation for the irreducible representaions Df{R) of the Sp group is 



1 



1 



This leads to the following relation for characters, true for any Q and Km. Sp] 

^ E xi{PQ)xARP-') = ^f.^xAQR)- (A6) 

The corresponding character relation for U (n) , 

/ idU)xfiUV)UWU-') = 5f,jXf{VW), (A7) 

is true for any V and W inU (n) . 

Given aU E U{n) and a P E Sp, define {U)p to be J^iUup, where the sum is taken 
over all the indices in the set / = (^1^2 ■ ■ - ip), each covering its full range of values from 1 to 
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n. The index set Ip as well as Unp are defined at the beginning of Sec. III. If P consists of 
ctj cycles of length i, then 

(U)P = n ■ (A8) 

j 

Since {U)p depends only on the cycle structure of P, it is the same for two P's in the same 
class. 

The crucial input to the computation of the integral is Frobenius formula, 

iU)p^j:XfiU)Xf{P)- (A9) 
/ 

Applying it to {UV)e and {WU~^)e, and using (A7) to integrate, one arrives at the expression 

/ {dU){UVUWU% - Yl fxfiVW), (AlO) 

J J (if 

where df = X/(e) has been used. Next, use (A6) to introduce the factor 

= 4 E Xf{P)x,{P-') (All) 
P- PeSp 

to the right-hand side of (AlO), and use the Frobenius formula again (note that Xg{P~^) — 
Xg{P))- This allows (AlO) to be written as 

[{dU){UVUWU^),= E j:^xf{P){vw)p 



Pes, f P^-df 
R,S'es, f {P^-^df 



E j:7^Xf{RS'){VW)RS'. (A12) 



Introducing the shorthand 

SkIj, SjL^, ^ ^kaiH(a) Kes'(a)^ ' (Al3) 

the sum over R and S' on the right-hand side is 

E xARS'WlkWij 5ki 5jl^,, = E xARS'WlkWij 5ki Sj^_,Ls, 

R,S' R,S' 

= E xARS'WlkWij^ 5ki Sjl,, 

R,S' 

- E XfiRS')VLKWi^_,j 6ki Sjl,, 

R,S' 

- E xARS'WlkWij 5ki^ 5jl,,. (A14) 

R,S' 
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Since V and W are arbitrary, the coefficients of VlkWu on both sides must be the same. 
Hence 

{I, J\K, L) = / {dU)UhUKL = T.ll t4iT SjL„Xf{RS'). (A15) 
•' f R,S' {P^-^df 

Let us now apply (A15) to the special case (/, J|/, Jq). Since K — I, the factor 5//^ 
is non-zero if and only if i? G Similarly, since L = Jq, the factor 5jLg/ equals to 
^JJqs'- Hence S — QS' must be in the invariant group Qj. Now the argument of in 
(A15) is RS' — RQ~^S. Since it! is summed over a group and so is S, we may replace 
R by R^^ and S by S~^. With this replacement, the argument of Xf is (SQR)"^, Using 
Xf{SQR) — Xf{{SQR)~^), we finally arrive at the formula shown in (Al). 

Formula (Al) can also be obtained from the following version of the Itzykson-Zuber 
formula 

f{dU) exp [pTt (m.UM^U^)] = Y: ^X/(Mi)x/(M2), (A16) 

where the sum is over all the Young's tableaux, and \ f\ — J2i fi the number of boxes in a 
particular tableau. Using the orthonormal relation (All), Frobenius's formula (A9) can be 
inverted to read 

/ {dU) exp [pTr {m,UM,U^)] = -^^T.Xf{R)Xf{S){M^U{M,)s. (A17) 
Next, use the formula 

Xf{R)XfiS) = ^^ExfiRPSP-'), (A18) 
I/I- p 

which can be derived from (A5), to combine the two characters into one. 



^TVJiKIT E Xf{RS'){MMM,)p-.s'P. (A19) 



[{dU)exp[pTr{M,UM,U^)] = ^ ^ Xf{RPSP-'){M,)R{M,)s 

/ df p,R,5 

p\f\dj 

Since {M2)s depends only on the class that S lies in, the last factor is thus independent of 
P. Therefore the sum over P yields only a factor |/|!. 
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Identifying terms proportional to (3^ on both sides, one gets 
Y: {I,J\K,L){M,)iK{M,)r.j^Y.T^T.Xf{RS')iM,)iK{M,Uj5Ki,Sj^^^ (A20) 

I,J,K,L f [PTdf R,S' 

which agrees with (A15). 
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